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DNA bending is biologically important for genome regulation and is relevant to a range of 
nanotechnological systems. Recent results suggest that sharp bending is much easier than implied 
by the widely-used worm-like chain model; many of these studies, however, remain controversial. 

We use a coarse-grained model, previously fitted to DNA’s basic thermodynamic and mechanical 
properties, to explore strongly bent systems. We find that as the end-to-end distance is decreased 
sufficiently short duplexes undergo a transition to a state in which the bending strain is localized at 
a flexible kink that involves disruption of base-pairing and stacking. This kinked state, which is not 
well-described by the worm-like chain model, allows the duplex to more easily be sharply bent. It 
is not completely flexible, however, due to constraints arising from the connectivity of both DNA 
backbones. We also perform a detailed comparison to recent experiments on a ’’molecular vice” 
that probes highly bent DNA. Close agreement between simulations and experiments strengthens 
the hypothesis that localised bending via kinking occurs in the molecular vice and causes enhanced 
flexibility of duplex DNA. Our calculations therefore suggests that the cost of kinking implied by 
this experiment is consistent with the known thermodynamic and mechanical properties of DNA. 


INTRODUCTION 

Strong bending of DNA plays a critical biological role 
in genome maintenance and regulation, most famously 
evident in DNA-protein complexes such as those formed 
in conjunction with the Lac repressor [1], the TATA bind¬ 
ing protein [2] and the nucleosome 0111. Experiments 
have confirmed that these large bending fluctuations are 
not an artefact of DNA-protein co-crystallization, and are 
observed in a variety of sensitive solution-phase assays 
[ 51111 . In the regime where DNA is not strongly bent, 
worm-like chain (WLC) models [TifTO] that treat DNA as 
a semi-flexible rod are sufficient to describe its configura¬ 
tional behaviour [HHIl]. It has been claimed, however, 
that strongly bent configurations can be reached much 
more easily than would be predicted by the WLC as para¬ 
meterized to small fluctuations. In particular, deviations 
from WLC behaviour have been reported in experiments 
on cyclization [HHis], multimerization HU, in DNA mini¬ 
circles [HUS], in a stressed ring system [20H23] and a 
molecular vice |24| . 

Much of the evidence for super-WLC flexibility, has, 
however, been contested. For example, in DNA cycliza¬ 
tion assays, cyclized molecules with ^ 100 base pairs (bp) 
[IllIS] were reported to be orders of magnitude more 
likely than predicted by the WLC model as normally 
parameterized. However, other authors have criticized 
the interpretation of the experimental results and have 
suggested that non-WLC behavior only begins at shorter 
lengths dSlESlESI- Similarly, apparent evidence of high 
flexibility of DNA on short length scales obtained via 
AFM [27] has since been attributed to artefacts of the 
curve-fitting algorithms |28| . 


The WLC model, which assumes a harmonic cost on 
the local degree of bending, favours the distribution of 
bending strain throughout a DNA molecule. However, 
above a certain degree of local bending the WLC model 
may break down because it becomes favourable for DNA 
systems to localize strain within small regions, allowing 
the remainder to relax and reducing the overall cost of 
adopting highly bent configurations. A disruption of the 
B-DNA structure (including broken base pairs and/or 
stacking interactions), or “kink”, is a candidate motif that 
would allow such enhanced flexibility and localization of 
bending [T8l[26l[29H32] . 

Nevertheless, although there is experimental evidence 
for kinking in strongly bent configurations (reviewed for 
example in reference [26]), it is not clear whether all of 
the proposed violations of WLC behaviour in the exper¬ 
iments listed above can be explained in this way. One 
reason for the lack of clarity is that a number of fun¬ 
damental questions about the thermodynamic cost and 
basic physics associated with the formation of kinks re¬ 
main open. It is therefore important to develop a sound 
theoretical underpinning of the underlying physics. Doing 
so would enable a systematic comparison between distinct 
systems and help establish whether experimental results 
are self-consistent and also consistent with basic DNA 
thermodynamics. 

Here we use oxDNA [33ti35] . a coarse-grained model of 
DNA, to address these questions. OxDNA incorporates a 
physical description of single- and double-stranded DNA, 
including the thermodynamics of hybridization and ba¬ 
sic mechanical properties such as persistence length and 
torsional modulus [33l|34]. As such, it is an ideal model 
for exploring strong bending, highlighting general prop- 
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erties and providing evidence as to whether the reported 
observations of enhanced DNA flexibility are consistent 
with kinking as implied by the known thermodynamic 
and mechanical properties of DNA. Further, as the model 
is generally applicable, oxDNA can be used to compare 
strong bending in a variety of distinct systems. Here, in 
the first of a series of papers investigating the biophysics 
of highly bent DNA, we consider a “molecular vice” m- 
In an accompanying paper we study cyclization [36], and 
in an upcoming future publication we consider minicircles 

1271 . 

In particular, in this paper, we first investigate the 
general thermodynamics of the strongly bent regime of 
DNA. Using this insight, we corroborate the results of 
recent “molecular vice” experiments [24|, which provide 
one of the most direct probes of strong DNA bending. 
Specifically, we reproduce the experimental high-salt buck¬ 
ling transition and approximate critical buckling length. 
We then directly probe the microscopic configurations re¬ 
sponsible for this transition, validating the experimental 
interpretation of a kink-induced buckling transition. 


MATERIALS & METHODS 
Concept of the “molecular vice” 

In an elegant set of experiments. Fields et al. deployed 
a “molecular vice” to probe the strong bending regime of 
DNA, providing evidence of a buckling transition at high 
salt that cannot be explained by a WLC model [24| . They 
infer kinking as the physical mechanism underlying this 
transition. The system they devised consists of a DNA 
hairpin and complement strand (Figure Q; the hairpin 
has two regions, a loop composed of A^oop nucleotides 
and a predominantly poly (AT) stem composed of 49 base 
pairs (39 A-T base pairs adjacent to the loop and 10 
terminal G-C base pairs). Complement strands of differ¬ 
ent lengths are hybridized to the hairpin loop to form a 
duplex of < Aioop base pairs. The hairpin acts as 
a force clamp, imposing a bending stress on the duplex. 
Formation of a duplex within the loop may require part 
of the hairpin stem to be opened, which requires work 
against a salt-dependent unzipping force. Stem unzipping 
is monitored via FRET (Fluorescence Resonance Energy 
Transfer), reporting the separation between the fluoro- 
phores attached to the first and last bases of the hairpin 
loop (Figure [^. 

This setup is powerful because it allows the observation 
of stem unzipping as a function of both the length of 
the complement strand and loop size. For example, at 
short there is minimal stress on the system - it is 

necessary neither for the duplex to bend nor the stem of 
the hairpin to unzip. At long (comparable to A^oop), 
the system cannot maintain full base-pairing without 
bending the duplex. At high salt Fields et al. observe 
for short loop lengths unzipping of the stem as A^j^p 
approaches A^oop base pairs. For some longer loop lengths. 


they then see a sharp re-zipping when A^j^p gets even 
closer to A^oop base pairs. Fields et al. argue that the 
high-salt behaviour cannot be explained by the WLC, 
and infer the formation of a sharp kink in the duplex that 
allows re-zipping |24) . 


The oxDNA model 

OxDNA [33H35] was developed to bridge the gap 
between all-atom molecular dynamics simulations and 
continuum or analytic approaches. It has been employed 
successfully for many systems, reviewed in reference [38] . 
Most importantly for this work, oxDNA has previously 
been shown to reproduce a range of structural transitions 
in response to applied stress [39H4T]. 

The philosophy of the model is to encode the well- 
understood and extensively characterized aspects of DNA 
physics, such as its basic thermodynamic, mechanical 
and structural properties, into a model with physically- 
motivated interactions. Briefly, each nucleotide (sugar, 
phosphate and base group) is modelled as a rigid body 
with three interaction sites. A complex set of anisotropic 
interactions between the nucleotides have been paramet¬ 
rized to reproduce structural, thermodynamic and mech¬ 
anical properties of single- and double-stranded DNA. The 
model is parameterized to a monovalent salt concentration 
of [Na^]=500 mM, where the electrostatic interactions are 
strongly screened. The thermodynamic, structural and 
dynamic properties of the model are fully characterized 
in references [33| and [34|. These include basic mechan¬ 
ical properties such as persistence length and torsional 
rigidity. 

The model has two parameterizations, average-base 
[33l [34] and sequence-dependent m In the average-base 
parameterization, the base pairing and stacking interac¬ 
tions of each base are identical; whereas in the sequence- 
dependent parameterization stacking and hydrogen¬ 
bonding interactions vary by base identity (e.g. A-T base 
pairs weaker than G-C base pairs). We use the average- 
base parameterization to elucidate the fundamental ther¬ 
modynamics of DNA bending since this illustrates the 
fundamental physics more clearly. Sequence-dependent 
thermodynamics are used for a more direct comparison 
with the molecular vice experiments and to elucidate 
sequence effects on the free-energy profile for bending. 


Bending simulations 

Simulations to probe the free-energy cost of bend¬ 
ing were performed with the virtual-move Monte Carlo 
(VMMC) algorithm (the variant in the appendix of Ref. 

[42] ) in the NVT ensemble at 298 K. As the entire DNA 
bending regime is of interest, including strongly bent 
states that are high in free energy, umbrella sampling 

[43] was employed to efficiently sample the free-energy 
landscape of DNA bending. The end-to-end distance 
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Figure 1. OxDNA representation of the “molecular vice” studied by Fields et al M- A duplex of A^up base pairs is formed by 
hybridizing the Aioop-nucleotide hairpin loop with a variable-length complementary strand. The duplex may be kinked (arrow) 
or unkinked; in the latter case there is an equilibrium between continuous bending in the duplex and unzipping in the stem. 
Fluorophores are attached to stem-end thymines (yellow) and are used for reporting hairpin separation via FRET measurements. 
For the representative system shown, the hairpin stem is 49 base pairs, A^oop = 36 and A^j^p = 30 bp. 


between the base interaction sites of the first and last 
bases along a strand was used as an order parameter 
for the umbrella sampling procedure, whereby the model 
Hamiltonian is augmented with an empirical potential 
that is iteratively adjusted to yield a roughly uniform 
sampling of the relevant values of R'^q. The statistical 
bias of this artificial potential is removed at the analysis 
stage. Further details of simulation procedure are given 
in Supplementary Section [SI A 


Molecular vice simulations 

Simulations of the molecular vice were conducted with 
the sequence-dependent parameterization of oxDNA us¬ 
ing the VMMC algorithm in the NVT ensemble at 25 °C. 
Error bars represent the standard error of the mean from 
6 independent simulations. These simulations were per¬ 
formed without any biasing through the use of umbrella 
potentials. Further details of simulation procedure are 
given in Supplementary Section jSl C 


To simplify interpretation, base pairing at the ends of 
the duplex were enforced through the umbrella potential 
to explicitly disallow “fraying” of duplexes. Fraying is a 
possible mechanism for relaxation of stress in bent systems 
[23l [24] , but we do not consider it in this section as we 
are interested in the bending of the main body of the 
duplex, not its ends. We allow fraying in the molecular 
vice system, but, as discussed by Fields et al [24], it is 
explicitly set up to inhibit fraying of the A^^^p base pairs. 
Results are essentially identical for any constraint in which 
1-5 nucleotides at either end of the duplex a re con strained 
to be base paired (Supplementary Section SIB); in this 
article we quote results in which we constrain only the 
final base pairs. Imposing such a constraint allows us 
to sensibly define and measure the end-to-end distance 
between the centre of the helical axis, Age, as the hrst 
and last base pairs are guaranteed to be present. Error 
bars in reported data represent the standard error of the 
mean from 5 independent simulations. 


Definitions 

Our study investigates the role of several structural 
defects, the detection and reporting of which rely on 
practical dehnitions. Namely, we require dehnitions of 
fraying, bubble formation and kinking. The identihca- 
tion of these defects requires, in turn, a dehnition of a 
formed base pair. Following previous work m, we dehne 
a base pair as formed if the hydrogen-bonding interaction 
between the two complementary bases is stronger than 
4.142 X 10“^^ J (^ 1 A;b ^ ^ = 298 K). We dehne fray¬ 

ing as the loss of base pairs at the duplex ends, while 
bubbles are dehned as the opening of base pairs in other¬ 
wise continuous stretches away from the duplex ends. 

A practical dehnition of a kink is less straightforward. 
Qualitatively, a kink is an area of pronounced curvature 
localized to a small region. There is some diversity of 
opinion on the quantitative description and classihcation 
of kinks [31]; however, they have in common the breaking 
of base pairs and/or stacking. Other requirements may 









4 


include a change in the relative orientation of nucleotides. 
Due to the ambiguity in the definition of a kink, we 
explore both structural and energetic definitions. 

Structurally, we define a kink using the relative orient¬ 
ation of consecutive nucleotides. In oxDNA, the orient¬ 
ation of each nucleotide is unequivocally determined by 
two orthogonal unit vectors, the base-backbone vector 
and the base-normal vector. The base-backbone vector 
connects the backbone and base interaction sites of each 
nucleotide. The model is designed such that, in a relaxed 
duplex, the base-normal vector at index i, a^, and at the 
consecutive index z -h 1 , are approximately parallel; 

that is, ~ 1 . A kink is defined to be present if 

• a^^i < 0 , a condition implying that the orientation 
of consecutive nucleotides differs by more than 90°. 

Energetically, we define a kink as a disruption in stack¬ 
ing between adjacent base pairs, where the stacking po¬ 
tential energy is less than one-tenth its typical value 
in a duplex (4.142 x 10“^^ J). For the oxDNA model, 
the structural and energetic criteria for a kink behave 
very similarly when kinks are reasonably probable due to 
strong bending of the duplex, as shown in Figure |^(c). In 
this regime, the detection of a kink is reasonably unam¬ 
biguous and straightforward as favourable kinks are ne¬ 
cessarily strongly bent, and both criteria do a reasonable 
job of identifying them. At larger end-to-end distances, 
the two criteria give distinct results, and caution should 
be used when interpreting the data; a fuller discussion is 
provided in Supplementary Section [SI D[ 




Ree j nm 


RESULTS ^ DISCUSSION 


DNA bending 


We start our discussion by characterising the basic 
thermodynamics of DNA bending in oxDNA. We simulate 
a 30-bp DNA duplex and report the free-energy profile 
AG(i?ee) as a function of the end-to-end distance i?ee 
(Figure [^(b)). The free energy of duplex bending has 
two distinct regimes, indicated by a change in slope at 
j^trans ^55 change in slope of the free-energy 

profile in Figure |^(b) correlates with a crossover from 
continuously bent to kinked macrostates (Figure He)). 
Representative oxDNA structures for the continuously 
bent and kinked 30-bp duplex (Figure |^(a)) reveal that 
duplex kinking typically induces a 1-3 bp bubble at the 
kink (further quantified in Supplementary Section S 2 A) 
and that this kink is loc ated n ear the centre of the duplex 
(Supplementary Section S 2 B). 

In Figure |^(b) we also compare oxDNA results to the 
WLC model; specifically to the analytic approximation 
to the WLC end-to-end distance developed by Becker et 
al. cni. This approach provides a good description of 
WLC behaviour across the near-rigid, semiflexible and 
flexible regimes. As expected, agreement between the 
WLC (when using the base pair rise and persistence length 
appropriate to oxDNA) and oxDNA is excellent for inter- 


Figure 2. (a) Representative OxDNA conhgurations for a 

30-bp duplex, depicting the kinked (arrow) and continuously 
bent states at Ree = 4nm and Ree = 8nm, respectively, (b) 
Free energy of bending AG{Ree) for oxDNA (solid black), with 
the decomposed free energies of the continuously bent (light 
green dashed) and kinked (light blue dashed) states overlaid. 
OxDNA is compared to the WLC model approximation of 
Becker et al. [TO] (red dotted). The comparison is fit-free, using 
known oxDNA parameters for persistence length (41.82 nm) 
and rise per base pair along the helical axis (0.34 nm) |34] 
Note that, rather than setting the minimum AG (Ree) to 
be zero as elsewhere in the text, we normalize AG (Ree) so 
that exp ( - ^^/k^T) = R(Ree) in both cases, (c) Probability 
of kinking Rkink using energetic (black) and structural (grey) 
criteria for a kink, which are based upon disruption of stacking 
and base orientation, respectively. 


mediate Ree- The deviations most pertinent to this paper 
arise at Ree ^ 6 nm where the oxDNA model develops a 
kink that allows the system to lower its free energy com¬ 
pared to the WLC model. Deviations also arise at longer 
Ree near the DNA contour length Rcontour, because the 
WLC neglects the finite intrinsic extensibility of DNA. 
Additionally, the WLC does not account for the finite 
thickness of DNA, which causes excluded volume effects 
at very short Ree, reported by oxDNA as a sharp rise in 
A G at Ree ~ 2 nm. 

According to the WLC model, the enthalpic cost of 
bending a duplex to some fraction of its contour length de- 
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Figure 3. (a) Length dependence of AG{Ree), the duplex 

DNA bending free energy, as a function of reduced end-to-end 
distance Ree /^contour • Error bars are comparable to the line 
thickness, (b) Length dependence of duplex DNA kinking. 
The structural criterion for kinking is shown; results using 
the energetic criterion are very similar in the kink transition 
region. 


creases with increasing duplex length. By contrast, there 
is no obvious reason why the free-energetic cost of forming 
a kink at the centre of the duplex should depend strongly 
on duplex length. We therefore expect that, as duplex 
length is increased, the DNA will need to be relatively 
more bent (smaller i^ee/^contour) to favour the formation 
of kinks. For sufficiently long DNA, we expect to see 
unkinked duplexes even for small i^ee/^contour- These 
expectations are confirmed in Figure long duplexes 
(> 73 bp) do not show any substantial kinking at any i^ee, 
and for shorter duplexes the length-dependence of the 
kink transition is as expected. The re-entrant behavior 
where Pkink drops again at extremely small i^ee/^contour 
arises from a change in geometry that is enforced by 
excluded volume effects; for a detailed discussion see Sup¬ 
plementary Section [82 C[ 

We next consider the question of what is the free-energy 
cost of kink formation, a question which turns out not 
to be so easy to answer unambiguously. One natural 
approach is to define a free-energy cost AGtrans that cor¬ 
responds to that at the middle of the kink transition, e.g. 
AGtrans = AG(i?gg^^®) = 18 A;b T for the 30-bp duplex 
in Figure This quantity is essentially a measure of 
work that must be done on the duplex by pulling its ends 
together before kink formation occurs. This numerical 
value is comparable to the estimate of Fields et al for 
“AGcollapse” EH, the bending stress released upon kink¬ 


ing for a similar system (in that kink formation is driven 
by pulling two ends of a duplex together). However, free 
energies must be carefully defined if quantitative compar¬ 
isons are to be made. In particular, neither AGtrans nor 
A Gcoiiapse is a generic free-energy cost of kink formation 
that can be directly compared with other systems in which 
kink transitions may occur. Furthermore, even when the 
bending stress is applied in the same way, AGtrans is not 
a constant, but, as we show in Supplementary Table [S^ 
for the current system, is a function of duplex length (be¬ 
cause the degree of bending at the kink transition depends 
on length). 

One might instead define a quantity AG^ink through 
the probability of spontaneous kink formation (in the 
absence of stress) at a specific location in the duplex. 
AGkink would in principle be directly measurable from 
kinking probabilities of relaxed duplexes. Although in 
principle this would be expected to be a generic (sequence- 
dependent) quantity, its value will be strongly dependent 
on t he crit eria used to identify kinks (Supplementary Sec¬ 
tion |ST^ . This is in sharp contrast to AGtrans, which 
is almost identical for both criteria used to identify kink¬ 
ing in this work. In the presence of sharp bending, the 
overwhelming majority of disruptions to the basic duplex 
structure are sharply bent regions with broken stacking 
interactions, which are unambiguously “kinks”. In an un¬ 
stressed system, however, the most common deviations 
from the basic duplex are not unambiguously identifiable 
as kinks, making the definition of a meaningful AGkink 
problematic. 

Even given well-defined criteria for identifying kinks, 
AGkink is not directly relatable to AGtrans- Two reasons 
for this difference are particularly evident in the context 
of bending duplexes. Firstly, to be a generic quantity, 
AGkink must be defined through the probability that a 
kink occurs that is centred on a specific base pair or (base 
pair step) in a duplex. In the context of duplex bending, 
however, kink formation can occur over a finite range 
of locations ( Supplementary Figure [S^, although kinks 
near to the centre are preferred, and this fact is refiected 
in AGtrans- Secondly, in pulling the ends of the duplex 
together to we do work on both the unkinked and 

the kinked state. Although kinks are more flexible than 
unkinked duplexes they still have a finite resistance to 
bending, due both to entropic effects and the need to 
strain the kink itself, as evidenced from the post-kinking 
slope in AG(i?ee) iu Figurej^b). To compensate, we must 
apply even more stress (than if the kink was completely 
flexible) to the system to reach the point at which kinks 
become favourable relative to unkinked duplexes, tending 
to increase AGtrans- 

That the bending stress at which kinking occurs in a 
specific system (AGtrans) is not equivalent to the intrinsic 
cost of kink formation (AGkink) emphasizes the difficulty 
in comparing between distinct experimental setups, in 
which the factors contributing to AGtrans niay well be 
different. For example, in some setups, it may be possible 
to relieve stress by kinking at one of a number of posi- 
















6 


tions, as observed in simulations of short minicircles m- 
In others, such as the molecular vice [24], the kinking 
location may be constrained to be near the mid-point of 
a duplex. The free-energetic cost of bending the kinked 
state may also be distinct in different systems. Unless 
care is taken, such differences will hamper quantitative 
comparison of free energies. Analytic models [TSlIM] used 
to extract free energies from data must therefore be care¬ 
fully chosen if quantitative comparisons are to be drawn. 
Indeed, one of the advantages of explicit simulation with 
coarse-grained models is that direct comparisons can be 
made between distinct systems at the level of physical ob¬ 
servables, rather than indirectly via free energies derived 
from analytic models. In this work, we will only make dir¬ 
ect comparisons of AGtrans for systems that are identical 
but for the DNA sequences, in which case differences will 
largely reflect differences in the intrinsic cost of kinking, 
^ ^kink • 

Given that the relative free energies of bending and 
kinking can also be perturbed with sequence variation 
|44| . we expect that sequence variation should affect the 
kink transition. We report sequence dependence in the 
free energy AG(i?ee) of bending for a 28-bp duplex (Fig¬ 
ure |4](a ,d)). It is clear from these data that the poly (AT) 
sequence undergoes kinking more easily than the poly(GC) 
sequence; AGtrans is lower, presumably due to the weaker 
base-pairing in the poly (AT) sequence. A secondary ob¬ 
servation is that the poly(AT) sequence is slightly more 
flexible prior to the kink transition in oxDNA. These two 
effects should, in principle, have opposing consequences 
for tending to increase and decrease 

spectively. It is clear from Figure [^(d), however, that 

the ease of kinking dominates and poly (AT) has a longer 

otrans 

-^ee 


We also observe that a 65 % GC sequence has a slightly 
lower AGtrans than the average-base model (Supplement¬ 
ary Table [S^. This behavior is expected since the duplex 
will kink and bend preferentially in AT-rich regions. Due 
to the lack of easily kinkable regions, the average-base 
model is harder to kink than would be naively expected 
from an ensemble of 50 % GC-content sequences. 

Having investigated sequence variation, we now turn 
our attention to a system in which we introduce a mis¬ 
match, a structural motif that substitutes a base pair 
with a non-Watson-Crick pair. It is reasonable to expect 
that kinking will localize to the mismatch, provided it 
is conveniently positioned, and AGtrans will be reduced. 
Following Fields et al [24], a mismatch was introduced 
at the mid-point of an otherwise intact duplex. We in¬ 
vestigated GG, GT, AA and CC mismatches. 

As expected, oxDNA reports that mismatched se¬ 
quences have a lower A Gtrans and undergo the kink trans¬ 
ition at longer (Figure [^(b,e)). The change in 

AGtrans Is consistent with typical free energy costs asso¬ 
ciated with mismatches [45]: oxDNA reports 16.1 A;b^ 
for the perfectly-matched duplex versus 12.8-13.5 A;b ^ fo^ 
the mismatches, a deviation of 3 h^T (see Table S2). 
Fields et al. also found evidence that ease of kinking 


was correlated with the sequence-dependent mismatch 
cost predicted by the SantaLucia model of DNA ther¬ 
modynamics [23]. OxDNA’s treatment of mismatches is 
extremely simple, neglecting any residual interactions or 
possible steric clashes between geometrically incompat¬ 
ible bases. Consequently, it does not strongly distinguish 
between different mismatches, and does not reproduce 
sequence-dependent mismatch variation. 

Like a mismatch, a nick, which is a break in an other¬ 
wise contiguous stretch of DNA whereby the backbone 
phosphodiester bond is severed, is a structural motif ex¬ 
pected to impact DNA bending. In the relaxed state, 
base pairing is preserved and coaxial stacking may occur 
between the adjacent, but now non-covalently bonded, 
nucleotides. It is relatively easy for a nicked system 
to give sharp kinks whilst only disrupting stacking and 
not hydrogen-bonding (Figure [^(f)). We observe, un¬ 
surprisingly, that the presence of a nick dramatically 
reduces the bending free energy at the kink transition 
(Figure 0(c)). In accordance with this lower bending free 
energy (AGtrans = ^-7 kQT), the crossover from continu¬ 
ously bent to kinked regimes is shifted to larger 
by about 2 nm. 

It is also worth noting that DNA with a nick is far 
easier to bend then normal duplex DNA after the kink 
transition: the post-kinking slope for the nicked system 
(Figure [^(c)) is much shallower than for the intact du¬ 
plex (Figure [^(b)). The absence of the constraint of the 
additional backbone bond in the nicked system makes 
extremely sharp bending very easy once coaxial stacking 
has been broken; this contrast highlights the work that 
must be done to strongly bend a kinked duplex without a 
nick. We note that for the nicked system, a more substan¬ 
tial difference between energetic and structural criteria 
for kinks can be seen (Figure [^(f)). This distinction 
arises because nicked systems undergo significant kinking 
at relatively large i?ee, when it is less advantageous for 
kinked sections to bend extremely sharply; hence, be¬ 
cause the structural criterion implies an effective angle 
> 90° between consecutive bases, it can underestimate 
the prevalence of kinks. 


Molecular Vice 

Having rigorously explored bending in a generic context, 
we turn to the “molecular vice” of Fields et al. [24], which 
provides an exquisite experimental probe of the behaviour 
of duplex DNA under bending stress. We use oxDNA 
to explore the microscopic configurations that underlie 
their experimental observations. We will show that our 
findings are compatible with the observations of Fields 
et al. at high salt concentrations. In our simulations, the 
thymines located at the first and last bases of the loop 
were chosen as the reporter for hairpin separation since 
they are the closest analogues to the experimental setup. 
Specifically, we define as the separation between 

these bases’ centres of mass. 
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Sequence effects 



Mismatch Nicked duplex 





Ree / nm 


Ree / nm 


Ree / nm 


Figure 4. (a) Sequence dependence of bending for a 28-bp intact DNA duplex, comparing the sequence-dependent and 

average-base oxDNA parametrizations. Error bars are comparable to the line thickness, (b) Impact of mismatches on the 
bending of a 28-bp DNA duplex, compared to intact DNA, using the oxDNA sequence-dependent parameterization. Mismatches, 
ordered by increasing free-energy cost, are as follows: GG, GT, AA and GG. Sequences from Fields et al. m- (c) Impact of a 
nick on the bending of a 28-bp duplex for oxDNA (solid black), with the decomposed free energies of the continuously bent 
(light green dashed) and kinked (light blue dashed) states overlaid, (d) Sequence-dependence of kinking using the structural 
criterion for a kink; results using the energetic criterion are very similar, (e) Impact of mismatches on kinking, (f) Impact of a 
nick on kinking, via energetic (black) and structural (grey) criterions. A representative oxDNA configuration of a nicked duplex 
with a kink (arrow) at Ree = 4nm is also depicted. 


We begin our discussion by considering the system with 
a loop length of = 36 nucleotides and a variable 

length duplex of A/(^up base pairs (Figure |^. OxDNA 
reveals distinct behaviour for small, medium and large 
values of A^^up- These behaviours are related to the most 
prevalent deformation modes from the canonical relaxed 
duplex DNA structure at each given length. Possible 
deformations to consider involve unzipping of the hairpin 
stem, and fraying, bending or kinking of the duplex. 

For short duplexes (A^^^p < 20) relative to the loop 

length Aioop = 36, the hairpin separation is 

comparable to that of a hairpin without an annealed 
complement strand (Figure [^(b)). Figure [^(c) and (d) 
show that kinking and fraying of the duplex, and hairpin 
unzipping, are small, and simulation snapshots such as 
Figure[^(a)(ii) do not show significantly bent duplexes. In 
this regime, the duplexes are too short to exert significant 
stress on the system. 

For intermediate-length duplexes (A^j^p = 20 — 28), 
^hairpin ij^creases due to stem unzipping (Figure [^(b)). 
The duplex is sufficiently long that the system cannot 
maintain full base-pairing without experiencing consid¬ 
erable mechanical stress. The system must either bend 
the duplex, or suffer a reduction in base pairing (either 
via unzipping the stem or fraying at the ends of the du¬ 
plex). In fact, the system adopts configurations in which 
the deformation is shared between continuous bending 
of the duplex and unzipping, as can be seen from Fig¬ 
ure [^(a) (hi). The bending energy released for each base 


pair unzipped decreases as the duplex gets straighter, 
whereas the cost of unzipping each base pair in the stem 
is approximately constant. The system can therefore min¬ 
imize its free energy by adopting a state which combines 
unzipping with moderate bending. Unzipping the hairpin 
stem is favoured over fraying the duplex for two reasons; 
firstly, stem unzipping contributes two nucleotides worth 
of single-stranded DNA to the hairpin loop (thereby re¬ 
ducing the tension), versus a single base for fraying the 
duplex. Secondly, the weak poly (AT) stem is easier to 
melt than the strongly binding GC-rich regions at both 
ends of the duplex (this was a deliberate design choice of 
Fields et al. [^). 

For long duplexes (28 < A^^^p ^ Ajoop), de¬ 

creases by re-zipping the stem and kinking, as reported 
in Figure ^(b-d); a typical kinked structure is shown in 
Figure [^(^(iv). Increasing A^^^p increases the stress in 
the system, leading to stronger unzipping and bending 
costs if the duplex does not kink. The cost of kinking, 
however, is much less strongly dependent on A^j^p (as 
was discussed in the general analysis of bending). Thus 
at sufficiently large A^^^p, it is favourable to form a kink 
in the centre of the duplex, relieving the bending stress 
elsewhere and reducing the need to unzip. 

As well as considering the behaviour at fixed A^oop, it 
is instructive to analyse the system as A^oop is varied, 
focusing on the highly stressed regime of A^^^p ^ A^oop 
(Figure [^(a) and (b)). We observe three regimes: For 
small Aioop < 30, unzipping and continuous bending 
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-^dup/bp 


Figure 5. Structural properties of the molecular vice as mod¬ 
elled by oxDNA for iVioop = 36 and variable iV^up- (sL)(i-iv) 
Representative structures for AT^up = 0? 10, 28, 34 bp, respect¬ 
ively. Stem-end thymines (yellow) are used for fluorophore 
attachment and reporting hairpin separation. A kink is high¬ 
lighted with an arrow, (b) is the separation of the 

two thymines bases at the interface between the hairpin stem 
and loop, (c) Pkink is the probability of kinking in the du¬ 
plex via the structural criterion, (d) Afray A^nzip are 
the number of disrupted base pairs via fraying in the duplex 
(green circles) and unzipping of the hairpin stem (dark green 
squares), respectively. 


are primarily observed, with almost no signal of kink¬ 
ing in ^hairpin gg ^ ^ gg^ 

we see a sharp 

kink transition in which kinked states become dominant 
^dup ^ ^loop (^loop = 30 sits near the crossover 
between these regimes). For larger A^oop, this transition 
is less pronounced, eventually fading away. Additional 
simulations at A^oop = 66 — 90, with A^jp^p = A^oop — 2, 
show the kinking probability decaying from ^ 30 % to 


—^ 20 bp — 26 bp —^ 30 bp —^ 36 bp 
—^ 40 bp —^ 46 bp —^ 50 bp —^ 56 bp 



Figure 6. Structural properties of the molecular vice for 
variable Aioop- (a) Kinking probability Pkink via the struc¬ 
tural criterion, (b) Hairpin separation Maxima in 

^hairpin highlighted (black circles), (c) The experimental 
FRET efficiency PpRET as reported in Ref. [24]. The available 
data is for Aioop = 30,36,40,46, 50. 


^ 2 % (Supplementary Figure [S7|). 

Why does a kink transition occur only for a finite range 
of Aioop? Unzipping of base pairs is extremely effective at 
relaxing stress for short duplexes: each unzipped base pair 
opened allows the duplex to straighten substantially. Thus 
at short lengths, stress is primarily relaxed by unzipping, 
duplex bending is small and kinking does not come into 
play. As the duplex gets longer, each unzipped base 
pair has less of a relaxing effect on the duplex, as it 
provides a smaller contribution relative to the overall 
length discrepancy between the intrinsic length of the 
duplex and the end-to-end distance enforced by the vice. 
Consequently, the bending of the duplex is enhanced, and 
eventually a kink transition is observed when this bending 
is sufficiently severe. For longer duplexes, the enthalpic 
cost of continuously bending to satisfy the separation 
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imposed by the vice drops with duplex length (this can 
be shown directly using the WLC expression of Becker 
et al. HOD, whereas kink formation does not get easier. 
Thus the driving force for kinking is reduced, and kinking 
is absent at sufficiently long loop lengths. 

Another notable point is that the maximum in 
prior to the kink transition (indicated by black 
dots in Figure [51(b)) drops with increasing A^ioop values 
of A^ioop which kinking is observable. Why should this 
be? As mentioned above, unzipping of stem base pairs 
is increasingly ineffective for larger systems, and bending 
of the duplex becomes increasingly easy. Therefore the 
degree to which the system unzips as Aj^up ^ ^loop is 


reduced for larger systems, and does not grow 

as significantly prior to kinking. 

We now compare the oxDNA results to those of Fields 
et al [24]. Unlike oxDNA, Fields et al could not dir¬ 
ectly observe kinked macrostates. Instead, they inferred 
kinking from the non-monotonic behavior of the FRET 
signal between dyes attached to the first two bases of the 
loop (a proxy for unzipping - Figure |^(c)). They defined 
a critical buckling length at high-salt, which they inter¬ 
preted as the onset of localized bending via kinking, as 
a local minimum in the FRET signal. We can define an 
analogous quantity in oxDNA, the maximum in 
for a given 

In their higher salt experiments. Fields et al saw evid¬ 
ence of unzipping and no kinking for = 30, identified 

non-monotonic behavior of FRET for = 36,40 and 

saw no substantial change in the FRET signal with Aq^p 
for larger values of A^QQp. Broadly speaking, these three 
regimes are consistent with our observations: kink trans¬ 
itions occur for intermediate values of A^oop, with pure 
unzipping at smaller values and no unzipping or kinking 
at larger values. Quantitatively, we observe the beginning 
of a kink transition for A^oop = 30, slightly earlier than 
Fields et al.^ and see evidence of kinking at larger values 
of Aioop. The Aq^p values at which kink transitions occur 
for Aioop = 36 and A^oop = 40 in oxDNA are within two 
base pairs of those inferred by Fields et a/., as can be 
seen in Figure |^(c) (clearly, however, inferring the ‘exact’ 
onset of kinking is not trivial from the experimental data). 

The quantitative consistency of results may be even 
stronger than immediately apparent. OxDNA is para¬ 
metrized at [Na^]=500mM, meaning that the stem of 
the vice should be somewhat more resistant to unzip¬ 
ping than in the experiments, which were performed at 
[Na^]=250mM, perhaps explaining the onset of kinking 
for Aioop = 30 that is not apparent in the work of Fields 
et al (f or add itional discussion, please see Supplementary 
Section S3B) Consistent with this argument. Fields et 
al. do not find evidence for kinking at all at much lower 
salt concentrations, when the stem of the vice is less 
stable. Slightly smaller values of A^^^p at which buckling 
occurs for a given A^qqp in oxDNA may also be explicable 
in this way. 

Our results for larger values of A^oop also suggest that 
although kinking occurs, it does so following an increas¬ 


ingly small amount of unzipping. This is actually con¬ 
sistent with the observations of Fields et al. — compare 
the drop in FRET prior to kinking for A^qqp = 36 and 
Aioop = 40. We posit that unzipping prior to kinking 
may be undetectable through the methodology of Fields 
et al. for A^oop ^ 40. There is a highly non-linear rela¬ 
tionship between FRET signal and ffuorophore separa¬ 
tion, which is roughly equivalent to ERET is 

most sensitive near the Forster radius, Rq = 5.7 — 6.7 nm 
for the Cy3B/Alexa647 ffuorophore pair used to probe 
stem unzipping. OxDNA reports < 3nm for 

Aioop ^ 46, suggesting that it is possible that kinking 
occurs in experiment at larger values of A^oop without 
sufficient change in FRET signal for it to be observable. 
OxDNA does suggest, however, that kinking is never ob¬ 
served when Aioop is large enough that A^j^p r\j Nloop 
duplexes undergo continuous bending. 


CONCLUSION 

In this paper we have provided a full characterization 
of the thermodynamics of DNA bending in a coarse¬ 
grained model, oxDNA. When weakly bent, oxDNA is 
well-described by the worm-like chain (WLC) model. 
However, if a sufficiently strong degree of bending is im¬ 
posed, oxDNA tends to adopt an alternative conformation 
in which the strain is localized in a small region of dis¬ 
rupted base-pairing and stacking interactions. The cost 
of such local kinks are compensated by the relaxation 
of the rest of the molecule. Such alternative modes of 
relaxation are not captured by the WLC model and lead 
to an enhanced probability of sharp bends. Although a 
kink is far more fiexible than canonical B-DNA, bending 
is not completely free. Conformations are more restricted 
than for a duplex with a nick, for example, which is not 
constrained by two continuous backbones. 

Our results enable us to identify the end-to-end con¬ 
straint that is sufficient to cause kinking (or localized 
bending). Short duplexes exhibit localized bending at 
relatively larger end-to-end distance than long duplexes, 
while sufficiently long duplexes (> 70 bp) do not undergo 
kinking even if the ends are held as close together as pos¬ 
sible; they are long enough to achieve such a conformation 
through reasonably gentle bending distributed through¬ 
out the molecule. Additionally, AT-rich sequences exhibit 
localized bending at longer end-to-end distance than GC- 
rich sequences, due to the lower cost of disrupting weaker 
interactions. 

Our analysis of kinking reveals that one must be ex¬ 
tremely careful in inferring free energies of kinking from 
specific systems. A natural definition of AG^ink would 
be through the probability of spontaneous kink forma¬ 
tion in relaxed DNA at a specific site - this should be 
a generic property depending on the local sequence, and 
would be the quantity measurable if relaxed duplexes 
were examined for kinking. If kinking is enforced through 
external bending of duplex DNA, however, the amount 
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of work that must be done on the duplex state to enforce 
a transition, AGtrans, is potentially quite distinct from 
AGkink- Reasons include the fact that work may also be 
done in bending the kinked state to reach the point of 
transition, and the fact that specific experimental systems 
(such as the molecular vice [24]) may favour kink forma¬ 
tion at certain locations whereas other systems (such as 
DNA minicircles m) may allow kinks anywhere. These 
considerations make direct comparisons of free-energies 
of kinking as inferred through simple analytic models 
difficult, hindering the comparison of kinking in distinct 
experimental systems. 

Computational models such as oxDNA have the vir¬ 
tue that they allow comparison between different experi¬ 
mental systems in which kinking has been hypothesized, 
as they can be directly simulated with the same model 
and experimentally comparable observables can be com¬ 
puted. In particular, physical factors such as the work 
done in bending a kink, or the preference for some kink 
locations over others in an experimental setup, should 
be naturally incorporated (at least at a semi-quantitative 
level). Coarse-grained models parametrized to reproduce 
basic DNA mechanics and thermodynamics therefore offer 
a valuable tool for establishing whether proposed obser¬ 
vations of kinking are consistent with each other and our 
current knowledge of DNA. 

To this end, we have simulated the “molecular vice” of 
Fields et al. [53|, a direct probe of DNA bending. In this 
experiment, the hybridization of a strand to a hairpin 
loop can lead to a range of behaviours interpreted by 
Fields et al. as duplex bending, kinking and hairpin stem 
unzipping dependent on the length of the strand and hair¬ 
pin loop. Our results are consistent with the observations 
of Fields et a/., support their interpretation that kinking 
is responsible for rezipping at high salt, and provide a 
more detailed mechanistic perspective on the data. We 
also propose that duplex kinking might be present for 
longer duplexes than can be directly inferred from the 
experimental data, due to the limits of FRET precision 
and minimal unzipping prior to kinking. 

OxDNA is a simplified model, and the level of agree¬ 
ment with Fields et al. should not be over-interpreted. 
Although reproduction of the data is evidence that the 
model is well-founded, a number of modelling simplifica¬ 


tions, and the slight difference in salt conditions in exper¬ 
iment and oxDNA, might give rise to small quantitative 
deviations. The important point is that the general agree¬ 
ment between the experimental data and oxDNA helps 
to further validate Fields et al. ’s inferral of kinking, and 
to show that this behaviour is consistent with the known 
thermodynamics and mechanics of DNA as captured by 
the oxDNA model. Furthermore, we have recently also 
performed careful comparisons between oxDNA simula¬ 
tions and experiments on DNA cyclization [36| and mini¬ 
circles m Reasonable consistency with experiment in 
these studies allows us to conclude that the experimental 
evidence of kinking in these systems is broadly consistent, 
and gives us further confidence that oxDNA describes 
the basic physics of kinking well. This consistency is 
particularly important given the controversy surround¬ 
ing the existence of kinks and enhanced DNA flexibility 
diiiH]. We note that our observation that kinking may 
go undetected for intermediate loop sizes in the setup 
of Fields et al. is potentially important in this regard; 
otherwise the absence of kinking for fairly short A^oop 
might be surprising when compared to digestion assays 
on minicircles HU IE], particularly as simple mechanical 
arguments suggest that the geometry of torsionally re¬ 
laxed minicircles is less conducive to kinking than the 
geometry of the molecular vice. 
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SI. EXTENDED METHODS 
A. Bending simulations 

Bending simulations were performed in three phases: 
exploratory, equilibration and production. In the explorat¬ 
ory phase, we iteratively adjusted the umbrella sampling 
biasing potential to yield a flat probability distribution 
P{Rqq)^ where is the end-to-end distance between 
the base interaction sites of the first and last bases on a 
single strand using a discrete potential with increments of 
0.4259 nm. As stated in the main text, the duplexes were 
constrained not to fray, allowing us to sample in terms of 
i?ee, the end-to-end distance between the centres of the 
helical axis at the first and last base pairs, choosing an 
Ree bin-width of 0.085 nm. Non-native base-pairing was 
also forbidden. As the iteration was performed manually, 
simulation times varied in the range 10^-10^ VMMC steps 
per particle. 

In the equilibration phase, we equilibrated the system 
for 10^ VMMC steps per particle using the aforementioned 
umbrella biasing potential. Adequate sampling (number 
of transitions between end-to-end distance bins) and de- 
correlation (via a block average decorrelation method) of 
potential energy, bubble size, fraying (where applicable) 
and structural kinking were checked. For the produc¬ 
tion phase, each production simulation (five independent 
simulations per measurement) was initialized with a stat¬ 
istically independent starting configuration and a unique 
random seed. This is accomplished by randomly draw¬ 
ing starting configurations from the equilibration phase, 
one for each production run, ignoring the first 2 x 10 ^ 
VMMC steps per particle of equilibration. The produc¬ 
tion run time for each simulation was 10 ^ VMMC steps 
per particle. For reference, the characteristic decorrela¬ 
tion time for the potential energy is r\j 10 ^ VMMC steps 
per particle, while the decorrelation time for kinking is 
^ 10^ VMMC steps per particle. The “seed moves” used 
to build clusters in the VMMC algorithm were: 

• Rotation of a nucleotide about its backbone site, 
with an axis chosen uniformly on the unit sphere, 
and with an angle drawn from a normal distribution 
with a mean of zero and a standard deviation of 
0.10 radians. 

• Translation of a nucleotide, where the displacement 
along each Cartesian axis is drawn from a normal 
distribution with a mean of zero and a standard 
deviation of 0.85 A. 

Additional simulations were performed to measure the 
free energy of kinked states for a 30-bp duplex at lar¬ 
ger values of i?ee- In this case, as well as biasing i^ee. 


the umbrella potential was used to forbid states without 
kinks (separate simulations were performed for each kink 
criterion). Umbrella potentials for Ree were obtained via 
manual iteration at a preliminary stage. Simulations were 
started from kinked configurations and thermalized for 
3.3 X 10^ VMMC steps per particle. Subsequent data 
collection lasted for 3.3 x 10^ VMMC per particle. Five 
independent simulations were performed. The VMMC 
“seed moves” differ slightly from those mentioned previ¬ 
ously: 

• Rotation of a nucleotide about its backbone site, 
with the axis chosen from a uniform random dis¬ 
tribution and the angle from a normal distribution 
with mean of zero and a standard deviation of 0.15 
radians; 

• Translation of a nucleotide with the direction chosen 
from a uniform random distribution and the dis¬ 
tance from a normal distribution with mean of zero 
and a standard deviation of 1.28 A. 


B. Constraints on fraying 


To isolate the effects of DNA bending, as opposed 
to alternative deformation modes such as fraying of the 
terminal base pairs, we impose a condition preventing 
fraying in our DNA bending simulations. Note that this 
constraint was not applied in simulations of the molecular 
vice system. Using an additional order parameter, the 
30-bp system was constrained such that the terminal base 
pairs remained hydrogen-bonded, thereby eliminating the 
possibility of fraying at duplex ends. We compared con¬ 
straints where 1, 3 and 5 base pairs at the duplex ends 
were prevented from opening, finding that the various 
conditions yield comparable results (Figure SI). 

To further probe the behaviour of the 30-bp system, 
in separate simulations we allowed the system to fray 
freely or preserved hydrogen bonding along the entire 
duplex. Unsurprisingly, as frayed base pairs can adopt a 
large range of conformations including many that reduce 
Ree, allowing fraying as an additional deformation mode 
reduces the bending free-energy penalty and delays the 
onset of kinking. Preserving hydrogen bonding through¬ 
out the duplex also delays kinking, but by limiting kinked 
states to those that do not involve disrupted base pair¬ 
ing; only those that break stacking but not base pairing 
remain, and thus the free energy cost of kinking is higher. 

Since a system in which fraying is permitted may not 
have terminal base pairs, the usual definition of Rqq re¬ 
porting the separation between the helical axis centres of 
the terminal base pairs, is not well defined. Therefore, for 
this system we instead use defined as the separation 
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(c) 




Ree or Rqq / nm 




Figure SI. The impact of constraining the terminal 1, 3 and 5 
base pairs from fraying on the free energy of bending AG{Ree) 
and kinking in an = 30 duplex. Limiting cases allowing 
fraying and enforcing hydrogen bonding along the entire duplex 
are also shown, (a) Free energy of bending. Probability of 
kinking based on the (b) energetic or (c) structural criteria 
for kinking. In all panels, for the frayed system Ree (the end- 
to-end distance between the centres of the helical axis at the 
first and last base pairs) is ill-defined, so the very similar Ree 
(the end-to-end distance between the base interaction sites of 
the first and last bases along the umbrella sampled strand) is 
used instead. 


between the base interaction sites of the first and last 
bases along a single strand. 


C. Molecular vice simulations 


A kin to the bending simulations of duplex DNA (Sec¬ 
tion SI A), the molecular vice simulations were performed 


in three phases: assembly, equilibration and production. 
An exploratory phase is not required since the simula¬ 
tions were unbiased, i.e. umbrella sampling is not used. 
As in the bending simulations, non-native base pairing 


was disallowed for simplicity. In the assembly phase, the 
system was initialized as two single strands. Harmonic 
traps were used to bring the stem ends, as well as the 
loop and complement strand, into proximity. Assembly 
with traps requires ^ 5 X 10^ VMMC steps per particle, 
significantly faster than free assembly. 

The assembled systems were equilibrated for 
10^ VMMC steps per particle. For the production 
phase, each production simulation (six independent 
simulations per point) was initialized with a statistically 
independent starting configuration and a unique random 
seed. This was accomplished by randomly drawing 
starting configurations from the equilibration phase, 
one for each production run, ignoring the hrst 2 x 10^ 
VMMC steps per particle of the equilibration phase. The 
length of each production run was r\j 10 ® VMMC steps 
per particle. 

Ah simulations were performed with the sequence- 
dependent oxDNA parameterization. Where available 
(Aioop = 30, 36,40,46, 50), the sequences are identical to 
those used in the experiments of Fields et al. [ST]; where 
not available , similar sequences are used. Ah sequences 
are written 5^ to 3^ 

The Fields et al stem sequence is used for ah 
systems: GCC CGG CGG CTT ATA AAA TTT ATT AAT TAT 
ATA TTT TAT TTA ATA TAA T-Loop. A complete list of 
loop sequences is available in Table [STj 


D. Operational definitions of kinking 


In the main text, we used two ‘kink detectors’ to identify 
the emergence of kinked states: one with an energetic 
criterion and one with a structural criterion. We saw that 
both measures located the transition to a kink-dominated 
ensemble at approximately the same Ree- To understand 
the properties of these detectors further, we performed ad¬ 
ditional simulations for a 30-bp duplex in which a second 
umbrella bias was used to forbid unkinked states (separ¬ 
ate simulations were performed for each kink criterion). 
In Figure S2, we report AG{Ree) (the free energy of the 
kinked stat^ versus Ree obtained from these simulations. 
Curves were normalized by manually matching to the 
data obtained for small Ree in the original simulations in 
which unkinked states were permitted - the zero of free 
energy is then the same as in ?? (a). 

Both detectors identify deviations from the B-DNA 
structure of the duplex that are extremely rare in the 
ensemble of typical states, as can be seen from Figure [S^ 
In snapshots of simulations around they also ap¬ 

pear to do a good job of distinguishing between clearly 
kinked and clearly unkinked structures. They are not, 
however, perfect reporters for conhgurations that would 
be intuitively interpreted as ‘kinks’. Firstly, there is al¬ 
ways a degree of arbitrariness involved in dehning a kink: 
is a disruption to the duplex structure of the type that 
allows enhanced bending a kink even if the net bend¬ 
ing angle around the disruption is less than 90°, or less 
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A/ioop Loop sequence 


20 

TAG 

CGA 

TAA 

GCT 

TGG 

TCA 

TT 























26 

TCC 

CAC 

CGA 

TAA 

GCT 

TGG 

TCA 

TGC 

CT 





















30^* 

TCG 

CCC 

ACC 

GAT 

AAG 

CTT 

GGT 

CAT 

GCC 

CGT 




















36^ 

TGC 

CCG 

CCC 

ACC 

GAT 

AAG 

CTT 

GGT 

CAT 

GCC 

CGC 

CGT 


















40^ 

TCC 

GCC 

CGC 

CCA 

CCG 

ATA 

AGC 

TTG 

GTC 

ATG 

CCC 

GCC 

GCCT 
















46^ 

TCC 

GCC 

GCC 

CGC 

CCA 

CCG 

ATA 

AGC 

TTG 

GTC 

ATG 

CCC 

GCC 

GCC 

GCC 

T 














50^ 

TGC 

CCG 

CCG 

CCC 

GCC 

CAC 

CGA 

TAA 

GCT 

TGG 

TCA 

TGC 

CCG 

CCG 

CCG 

CCC 

GT 













56 

TGC 

CGC 

CCG 

CCG 

CCC 

GCC 

CAC 

CGA 

TAA 

GCT 

TGG 

TCA 

TGC 

CCG 

CCG 

CCG 

CCC 

GCC 

GT 











60 

TCC 

GCC 

GCC 

CGC 

CGC 

CCG 

CCC 

ACC 

GAT 

AAG 

CTT 

GGT 

CAT 

GCC 

CGC 

CGC 

CGC 

CCG 

CCG 

CCT 










66 

TCG 

CCC 

GCC 

GCC 

CGC 

CGC 

CCG 

CCC 

ACC 

GAT 

AAG 

CTT 

GGT 

CAT 

GCC 

CGC 

CGC 

CGC 

CCG 

CCG 

CCC 

GCT 








70 

TGC 

CGC 

CCG 

CCG 

CCC 

GCC 

GCC 

CGC 

CCA 

CCG 

ATA 

AGC 

TTG 

GTC 

ATG 

CCC 

GCC 

GCC 

GCC 

CGC 

CGC 

CCG 

CCG 

T 






76 

TCC 

GGC 

CGC 

CCG 

CCG 

CCC 

GCC 

GCC 

CGC 

CCA 

CCG 

ATA 

AGC 

TTG 

GTC 

ATG 

CCC 

GCC 

GCC 

GCC 

CGC 

CGC 

CCG 

CCG 

GCC 

T 




80 

TGC 

CCG 

GCC 

GCC 

CGC 

CGC 

CCG 

CCG 

CCC 

GCC 

CAC 

CGA 

TAA 

GCT 

TGG 

TCA 

TGC 

CCG 

CCG 

CCG 

CCC 

GCC 

GCC 

CGC 

CGG 

CCC 

GT 



86 

TCG 

CGC 

CCG 

GCC 

GCC 

CGC 

CGC 

CCG 

CCG 

CCC 

GCC 

CAC 

CGA 

TAA 

GCT 

TGG 

TCA 

TGC 

CCG 

CCG 

CCG 

CCC 

GCC 

GCC 

CGC 

CGG 

CCC 

GCG 

CT 

90 

TGC 

CGC 

GCC 

CGG 

CCG 

CCC 

GCC 

GCC 

CGC 

CGC 

CCG 

CCC 

ACC 

GAT 

AAG 

CTT 

GGT 

CAT 

GCC 

CGC 

CGC 

CGC 

CCG 

CCG 

CCC 

GCC 

GGC 

CCG 

CGC CGT 


Table SI. Loop sequences for the molecular vice system. Terminal loop thymines, used to compute are highlighted 

(red). 

Mismatch locations as used by Fields et al. |S1| are highlighted (blue). Following the experimental design, the mismatch is 
introduced by changing the sequence of the complement strand. 

^ Sequences used by Fields et al m- 


than 45°? Of our two detectors, the structural criterion 
is more exacting in this regard, which explains why it 
reports slightly fewer kinks at i?ee ~ than the 

energetic criterion. Although the two detectors do not 
agree exactly, however, they are reasonably consistent 
at Ree ^ is because duplex disruptions are 

only common at Rqq < because they allow large 

bending. We do not observe a significant number of con¬ 
figurations that meet the energetic criterion without large, 
localized bending because such states do not help to relax 
the stress applied, and hence remain rare. Thus most 
configurations that meet the energetic criterion meet the 
structural criterion (although a few are not quite bent 
enough). At i?ee > number of structures 

that satisfy the energetic criterion but not the structural 
criterion is larger, partly because less bending must be 
localized at kinks to relax the system elsewhwere. Note 
that for the nicked system, kinking occurs at much larger 
Ree due to the lower cost; hence systems can relax the 
bending stress without meeting the structural criterion, 
leading to a bigger discrepancy between the two detectors 
at Ree ~ 

There is a second reason why the curves for the ener¬ 
getic and structural detectors diverge at Ree > 
which is related to a second reason why the detectors are 
not perfect reporters for kinking. As can be seen in Fig¬ 
ure the free energy curves for systems constrained by 
the energetic and structural detectors exhibit a noticeable 
change in slope as Ree gets larger. This change in slope 
is analogous to the change in slope of the free energy 


of the unconstrained system, also shown in Figure |S2[ 
Such a slope change is characteristic of a crossover from 
one typical conformation (with a certain fiexibility) to a 
qualitatively different one with a distinct fiexibility. For 
our unconstrained system, the slope change is simply the 
transition between kinked and unkinked states. Thus at 
Ree ^ 10 nm, the typical states that satisfy the energetic 
and structural criteria for kinking are qualitatively dis¬ 
tinct from those that are common at Ree ^ 
example, the energetic criterion can be met by stretching 
a base pair step, and the structural criterion can be met 
by flipping out a single base from the helix and rotating 
it; configurations of this kind are more probable than 
those that we would intuitively describe as kinked if the 
duplex is not sharply bent These configurations highlight 
an additional difficulty in unambiguously defining a met¬ 
ric and cost for kinking. However, such states are only 
common relative to other configurations that satisfy the 
kink criteria at large Ree] they are still rare compared to 
normal states at large Ree and remain rare at all values 
of Ree as they do not help to alleviate bending stress. 


In summary, despite the ambiguities in defining and 
identifying kinks, only duplex disruptions associated with 
localized bending through large angles are common at 
Ree ^ intact duplex. Thus the two kink de¬ 

tectors are sufficiently consistent and provide meaningful 
results. 
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Figure S2. Free energies as a function of Ree for a 30-bp 
duplex constrained to kink by either the energetic or structural 
criteria. The free energy is normalized by manually matching 
to the data reported in the main text (Figure [^(a)) for kinked 
states at small Ree- Error bars are comparable to line thick¬ 
ness. Also shown for comparison is the data for a system not 
constrained to show kinking. 



Figure S3. Normalized histogram for the probability distribu¬ 
tion (^bubble) of fhe total bubble size in base pairs (Abubble) 
given the presence of a kink via either the energetic (black) 
or structural (grey) criteria for A^bp = 30 at (Ree) ~ 5.3 nm 
from the umbrella potential bin Rqq = 4.26 — 4.69 nm. Error 
bars are of the order of the point size. 


(a) Type-I like 


(b) Type-II like 






Figure S4. OxDNA representations of type-I like and type-II 
like kinks in 30-bp duplex. Two perpendicular views are given 
in each case, (a) Type-I like kinks have disrupted stacking 
with preserved base pairing, (b) Type-II like kinks disrupt 
both stacking and base pairing. 


tion of both base pairing and stacking). Due to its coarse¬ 
grained nature, we have less confidence that oxDNA has 
sufficient molecular detail to contribute reliably to this 
discussion; furthermore, the distinction between type-I 
and type-II kinks are of less direct relevance to this study. 
We therefore do not focus on the classification of kink 
subtypes, although we note that we see config ura tions 
which could be assigned to both types (Figure [S4| ) and 
that most kinks involve broken base pairs (i.e. are of 
type-II) - for example, only 10% of the kinks for the 
states contributing to Figure [S3] only involve disruption 
of stacking. 


B. Localization of kinks in bending simulations 


S2. ADDITIONAL ANALYSIS OF DNA 
BENDING 

A. Bubble size distribution 

The distribution of bubble sizes reveals the frequent 
presence of a 1-3 bp bubble in kinked duplexes (FigurejS^. 
Briefly, bubbles are defined as loss of base pairs away from 
the duplex ends, while loss of base pairs at duplex ends is 
classified as fraying. Unsurprisingly, bubbles are almost 
always co-localized with kinks, so we take total bubble size 
along the duplex as a proxy for bubbles at a kink. Results 
for total bubble size given the energetic and structural 
criteria for a kink are very similar. 

Nomenclature has been developed in all-atom MD sim¬ 
ulations |S2HS4| to distinguish between type-I (disruption 
of stacking but not base pairing) and type-II kinks (disrup- 


When the ends of a duplex are pulled together, it would 
be expected that kinking occurs near the midpoint of the 
duplex. Kinks far from the midpoint would result in 
incommensurate relaxed duplex sections, which are then 
less likely to have their ends in close proximity. Analysis 
of the kink probability distribution confirms that kinks 
are localized near the midpoint (Figure [S^. 

Since structural and energetic kinks are a property 
defined between pairs of consecutive bases (see Defini¬ 
tions section of main text), we report the kink location 
distribution in base steps, where base step i refers to a 
kink between base index i and z + 1. For every sampled 
state, all base steps along one strand that satisfy the 
relevant kink criterion are included in the histogram. To 
normalise, we divide by the total number of detected 
kinks along both strands in all sampled states. Kinks are 
considered distinct if they are separated by more than 6 
base steps. Thus the normalized kink distribution gives 
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the probability that a given pair of neighbours contrib¬ 
utes to an arbitrarily chosen kink. To compare multiple 
lengths, we report the probability of kinking versus the 
adjusted base pair index, centred around the mid-point 
of the duplex, {N\^p — 2)/2. The —2 offset accounts for 
zero-based numbering and indexing the base step to base 
index i instead of z + 1. 



base pair index (adjusted) 


Figure S5. Normalized histogram for the probability distri¬ 
bution of kink localization for the energetic and structural 
criteria for intact (ATp = 20, 30,41, 52) and nicked (ATp = 30) 
duplexes. We sample all duplexes at constant end-to-end dis¬ 
tance from the umbrella potential bin Rqq = 4.26 — 4.69 nm. 


For the intact duplexes, the distribution is quite wide, 
^ 10 base pairs, suggesting that entropic effects are quite 
prevalent and some degree of incommensurately can be tol¬ 
erated. The distribution also appears to become slightly 
broader for longer duplexes. Importantly, the preference 
for localisation about the midpoint of the duplex is be¬ 
cause this is the position that provides the maximum 
decrease in end-to-end distance: in relaxed duplex DNA, 
or in DNA minicircles, one would not necessarily expect 
a preferred site for kinking with the average-base para- 
meterisation of oxDNA. 

For the nicked duplex at A^p = 30, the distribution is 
quite narrow, ^ 2 base pairs, reflecting the strong pref¬ 
erence for localisation of kinking at the nick (Figure [S^. 
Contrary to the intact duplexes, the preference for kink¬ 
ing at a nick is not just due to the bending stress imposed 
on the system: nicked sites have a lower free energy cost 
for kinking because base pairing need not be disrupted 
for a kink to form; only stacking must be compromised. 


(a) -Ree/-^contour ~ 0.1 


(b) R ee / -^contour ~ 0 




Figure S6. Representative OxDNA configurations for a 52-bp 
duplex at small Ree/Rcontour- (a) A kinked state, with Ree 
small but greater than 2 nm. (b) An unkinked, minicircle-like 
state with Ree ^ 2nm. With this constraint, a single kink 
cannot relax the majority of the bending stress, unlike in case 
(a). 


C. Re-entrant continuous bending 

At i^ee/^contour ~ 0, kinking appears to be less effect¬ 
ive at relieving bending stress; for all but the shortest du¬ 
plexes, there is re-entrant behaviour from kinking back to 
continuous bending. A typical kinked structure is shown 
in Figure [^(a); it consists of two long (relatively relaxed) 
duplex sections and a highly bent kink region. Given that 
Ree is measured in terms of the distance between terminal 
base pair centres, however, such a configuration cannot 
bring the ends into extremely close proximity (Ree < 2 nm, 
the approximate duplex diameter) due to steric effects. 
Instead, the blunt ends of the duplex must face each other 
to achieve Ree < 2 nm and this fundamentally changes the 
constraints placed on the rest of the system, and makes 
kinking opposite the blunt ends less effective at relieving 
the stress. 

The DNA duplex adopts minicircle-like configurations 
for extremely small values of Ree, as shown in Fig- 
ure[^(b). Minicircles can, of course, undergo kinking; 
however, the geometrical constraints mean that it is im¬ 
possible to localise the majority of bending stress into 
a single kink |S5| . unlike in systems with a slightly lar¬ 
ger separation of the blunt ends. Thus we obtain the 
seemingly paradoxical result that kinking is more likely 
for systems that have some flexibility in their ends than 
for systems, such as minicircles, where the ends are con¬ 
strained to be in extremely close proximity. 


D. Kink transition free energy 

For clarity, we present the bending free energy and end- 
to-end separation at the mid-point of the kink transition 
for all simulations in Table |S2j 
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A frtrans 

/ keT 

otrans / 

Aee / nm 

Length / bp 

20 

15.9 

4.7 

30 

18.1 

5.5 

o 

CO 

6.0 

8.0 

41 

19.7 

5.3 

52 

21.4 

4.3 

Sequence / 28 bp 

poly(GC) 

21.1 

5.1 

Avg-base 

18.0 

5.4 

poly(AT) 

14.8 

5.8 

Mismatch / 28 bp 

Intact 

16.1 

5.7 


13.0 

6.0 

mGT 

13.3 

6.0 

“AA 

12.8 

6.1 

“GG 

13.5 

5.9 

Nicked / 28 bp 

“poly(GC) 

7.9 

7.2 

^Avg-base 

6.7 

7.3 

"poly(AT) 

4.8 

7.6 


Table S2. The bending free-energy AGtrans and the end- 
to-end distance i^trans at the midpoint of the kink transition, 
which we define as the point at which kinks are present in 
50 % of states according to the energetic detector, for all the 
different systems studied here. Results for 62 and 73 bp are 
not given because although the kinked microstates occur, the 
probability of kinking never reaches 0.5. 

^ Nick at the mid-point of the duplex. 

™ Single mismatch in the duplex with the same sequences as 
used by Fields et al. m 


S3. ADDITIONAL DATA FOR THE 
MOLECULAR VICE 

A. Long-length limit 

In the main text, we show the structural properties 
for variable loop length {N\qq^) and duplex length (A/^^up) 
(??). The data suggest that kinked configurations are 
gradually superseded by continuously bent states as A^ioop 
increases. This behaviour is also observed in our investig¬ 
ation of length-dependence in the free-energy of bending 
of intact duplexes (??), and would therefore be expected 
for the molecular vice. 



A^dup/bp 


Figure S7. Length-dependence of (a) hairpin separation 
^hairpin, probability of kinking Pkink and (c) fraying 
(Affray) and stem unzipping (Aunzip) for the molecular vice in 
the Adup = (Aioop - 2) limit. 


We make the length-dependence of kinking, hairpin 
separation, fraying and stem unzipping in the molecular 
vice more explicit in Figure [S7| by profiling the longest 
Adup at a range of A^oop lengths. At the lowest values 
of Aioop 5 unzipping (and a small amount of fraying) is 
the dominant method of relieving stress. At intermediate 
Aioop, kinked states dominate. A transition from kinked 
to continuously bent states is observed at long Ad^p (and 
therefore long Aioop), with the midpoint at Ad^p = 58 — 
64 bp. By Adup = 88bp, duplexes are almost exclusively 
continuously bent (Pkink < 2%). Unsurprisingly, the 
length associated with this crossover is consistent with 
our bending simulations (??). 


B. Temperature-dependence 

We cannot currently study the exact salt concentration 
used by Fields et al. |S1| with oxDNA. However, we can 
probe the impact of increasing temperature, which may 















7 



•^dup/bp 


Figure S8. Temperature-dependent structural properties of 
the molecular vice for = 30 at T = 25 °C (black circle) 

and T = 45°C (gray square). 


also be somewhat analogous to decreasing salt concentra¬ 
tion, as both serve to weaken base pairing; in particular, 
reducing the unzipping force associated with the hair¬ 
pin |S6| . We also note that temperature is not perfectly 
controlled in gel-based experiments due to electrophoresis- 
induced heating of the gel running buffer. 


Temperature-dependence was explored with additional 
sim ulat ions at 35 °C (data not shown) and 45 °C (Fig¬ 
ure S8). We find negligible kinking at T = 45 °C for 


A/ioop = 30, consistent with the reasoning that weakened 
base pairing in the stem leads to easier unzipping. 
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